Lattice analogy of area-difference elasticity model for lipid-detergent bilayer 

vesiculation 
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The vesiculation process was examined in the hpid-detergent solution (dimyristoyl-sn-glycero- 
phoshatidylcholine/octaethylleneglycol n-dodecyl ether/water), using small-angle neutron scatter- 
ing experiments Q. When observing vesiculation proceeds from rod- like micelles to unilamellar 
vesicles, the transformation is induced by jump-like temperature increase and a monotonic increase 
in detergent concentration. Our numerical estimations of the vesicle shape parameters (the elastic- 
ity coefficients and its fraction on a macroscopic scale) are based upon the area-difference elasticity 
model P,!3|. Thus, we composed the numerical Monte Carlo method, which connects the macro- 
scopic and microscopic scales by the concept of self-avoiding random surfaces. 

PACS numbers: 02.70.Uu, 05.50.-|-q, 61.30.St, 61.30.Hn, 61.30.Cz 



INTRODUCTION 

We study the lyotropic processes, which are widely 
distributed in biological systems and define the self- 
organization phenomena. In lyotropic hquid crystals, the 
numerous examples (or even, in its clear majority), when 
topological genus of mesophase surfaces was changed at 
the phase transitions, are known. 

As Kiselev observed in small-angle neutron scat- 
tering investigations [l|, self-assembled vesiculation 
occurs between rod-like micelles and unilamellar vesi- 
cles with closed surfaces. Self-assembled vesiculation 
is induced by temperature jump of 100 K ■ min^^ 
from room temperature and increases as function of 
detergent concentration. The typical bilayer-vesicle 
shape dimyristoyl-sn-glycero-phoshatidylcholine/water 
{DMPC/H2O) was observed experimentally in ^. 
The different types of mesomorphism are distinguished 
by the nonionic surfactant influence in dimyristoyl- 
sn-glycero-phoshatidylcholine/octaethyllene-glycol 
n-dodecyl ether/water {DMPC/C12ES/H2O) and the 
drastic temperature jump. 

In principle, certain closed stationary shapes (pears, 
prolates, oblates, stomatocytes and nonaxisymmetric 
shapes with the symmetry D2h) are possible |2]- The- 
oretical models, which combine phenomenological and 
microcanonical methods were developed to interpret the 
spontaneous vesiculation phenomena in mesomorphism. 
Based on Helfrich theory, a spontaneous curvature cq 
(shape independent cq with discontinuous budding tran- 
sition) model and area difference model (shape indepen- 
dent A^ with continuous one). The area-difference elas- 
ticity (ADE) model 3] describes the energy as a sum of 
two terms: one depends on curvature at any point on the 
surface, and second depends on curvature-induced area 
difference between the inner and outer monolayers. In 
the bilayer-couple model by Svetina et al. [S| , the bilayer 
structure is modelled by representing its two monolayers 



as closed neutral surfaces with a constant separation dis- 
tance. The elastic energy is composed of the local and 
nonlocal bending energies of the membrane. 



MODEL OF SHAPE TRANSFORMATIONS 

At a spasmodic rise in temperature, a rod-like micelle 
becomes unstable and transforms into a bilayer, which 
then collapses into a hollow sphere. Disintegration of the 
rod-like micelles and the bilayer were not experimentally 
observed [ij, presumably, due to a fast jump. All of these 
phenomena look like a first order phase transitions, and 
there are situation when the symmetries of both phases 
are not in some ratio. 

The typical budding process has been phenomenologi- 
cally described &\ and the shape equation for the axisym- 
metric equilibrium shape was obtained from some general 
variational ansatz. In a pure DM PC solution, generative 
only the lamellar phase, it is found the lateral compres- 
sion modulus and bilayer bending rigidity. However, at 
addition of nonionic surfactant C'i2Es, it becomes diffi- 
cult to analytically predict the same parameters of the 
new mesomorphic series. Using the ADE model pltzl it is 
possible to determine the minimum energy shape of the 



vesicle as a function of the reduced volume v 



, where 



V is the enclosed volume and Vq is the initial volume. 
The reduced volume v is concerned with the measured 
hvdrodynamical radius of a vesicle by certain relations 
[2|. To define the vesicle structure parameters, methods 
that combine the Helfrich spontaneous curvature theory 
0,13 and the bilayer coupHng model were developed. 

In terms of a microscopic canonical ensemble and lat- 
tice Ising-like models, the main mesophases structure pa- 
rameters may be connected by a spontaneous curvature, 
which is constructed on dual lattices jSj by self-avoiding 
random surfaces. One would think so, that in frame of 
any spin variables, the mixture content is disregarding. 



LATTICE MONTE CARLO METHOD 

As noted by Caselle |8j, a gas of self-avoiding ran- 
dom surfaces (SARS) with unconstrained topology and is 
weighted by the usual area term e~^^^ where yl is a sur- 
face area and (3 = l/fc^T (fc^ is the Boltzmann constant, 
T is the absolute temperature), belongs to the same class 
of the Ising model. As other kinds of coupling are intro- 
duced, it is possible to generate these surfaces on the 
body centered cubic {bcc) lattice and the dual lattice, 
like the extrinsic curvature. The reduced Hamiltonian of 
the Ising model has the form 



/3iJ: 
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with the partition function Z(f3s, Ph, Pt) — J2a=±i ' 
where /3s, Ph, Pt are the coefficients of couplings, < ij > 
and << kl >>, (ijkl) are three kinds of coupling, and 
(7i=-l,0,l are the spin components, which correspond to 
the different types of mixture molecules (lipid, surfactant 
and water respectively), ay — <Ji<jjak<Ji- The coupling 
coefficients are normalized by action A in terms of the 
Riemannian curvature R i, which is in accordance to 
the Gauss-Bonnet theorem J (P^^/gR = ^^ Ry = x{S)- 
Ri is the curvature components for four site spin config- 
urations on the hcc lattice. The concentration term Pt 
in this model is in respect to by the chemical potential fx 
[9| . This treatment allow the numerical estimations to be 
conducted with the lattice Monte Carlo technique for the 
free energy and structural characteristics of the vesicular 
phases. 

According to these calculations, all of transitions with 
qualitative surface changing genus, at Pt range between 
0.4 and 0.8, may correspond to the appropriate transi- 
tions, found in [g] (Fig. 1.). The structure parameter of 
a vesicle shown at Fig. 2. 

The shape function r{9, </?) corresponds to the surface 
radius Rg for the other follows numerical simulations. 
This model conveniently calculates the including volume 
and other main structure parameters to plot the phase 
diagrams of lyotropic mixtures. 



SUMMARY 

The ADE model satisfactorily describes the vesicle 
phase transformations, where the non-homogeneous lo- 
cal properties can be expressed by the elasticity mod- 
ulus and its fraction in a macroscopic scale. However, 
the presented model is well suited for studying numer- 
ous biophysical phenomena and complements the exist- 
ing conceptions. In addition to phenomenological ADE 
calculations, the lattice treatment of the lyotropic phase 
behavior of the aggregates is essential for directly control- 
ling the concentration, which is a determinative param- 



lattice tern perature in terms o1[9] 



FIG. 1: The heat capacity dependence of temperature pa- 
rameter P for binary mixture with nonionic surfactant, Monte 
Carlo error is 10"'* 
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FIG. 2: Monte Carlo simulations for vesicle structural param- 
eter (mean curvature R(A)) evolution at 303 K 



eter in the lyomesomorphic transformations. Computer 
cluster simulation analysis is useful to investigate the lo- 
cal curvature properties of different aggregates, including 
vesicles. 

The author thanks to M.A. Kiselev for explaning of 
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O.P. Santillan and the anonymous referee for useful cor- 
rection of English style. 
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